Solving the Observer Metric 
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The analysis of modern cosmological data is becoming an increasingly important task as the 
amount of data multiplies. An important goal is to extract geometric information, i.e. the metric of 
the cosmos, from observational data. The observer metric is adapted to the reality of observations: 
information received along the past null cone, and matter flowing along timelike lines. It provides 
a potentially very good candidate for developing a general numerical data reduction program. As 
a basis for this, we elucidate the spherically symmetric solution, for which there is to date single 
presentation that is complete and correct. With future numerical implementation in mind, we give 
a clear presentation of the mathematical solution in terms of 4 arbitrary functions, the solution 
algorithm given observational data on the past null cone, and we argue that the evolution from one 
null cone to the next necessarily involves integrating down each null cone. 
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I. INTRODUCTION 

Cosmology is all about understanding the observed 
universe, and since Einstein's field equations are central 
to that endevour, the primary problem is to determine 
the cosmic geometry, i.e. the metric, from observations of 
its matter content. A full understanding of the dynamics 
of the universe cannot be separated from understanding 
its geometry. Historically, it was very difficult to deter- 
mine cosmological data with any precision, and the as- 
sumption of a homogenous universe, which allowed a sim- 
ple metric form, was entirely sufficient and indeed very 
fruitful. Consequently the problem was reduced to one of 
finding the best-fit parameter set, rather than determin- 
ing the metric. Recent decades have seen a considerable 
improvement in the quality and quantity of cosmological 
data, mapping much more accurately the matter distri- 
bution and the structures that exist. Therefore, relaxing 
the assumption of homogeneity has become an important 
task. 

The idea of determining the spacetime metric from 
observational data was first investigated by Kristian 
and Sachs f^, and followed up by Ellis, Stoeger and 
others in an important series of papers that consti- 
tute the "observational cosmology" (OC) programme 
0, i i, i, i 0, 1 i E, Lll, 12, ij, tii. In this pro- 
gramme, they introduced observer coordinates based on 
the past null cones of a single observer's worldline, an 
idea originally due to Temple [lB|, because traditional 
time and space coordinates are not well adapted to cos- 
mological observations. They also introduced the 'fluid- 
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ray tetrad' |3|, including a set of spin coefficients, and 
their basic equations are derived from this formalism. 
Although a general form has been given for the observer 
metric, work has concentrated on the spherically sym- 
metric case, and to a lesser extent its perturbations. 

There has also been parallel work relating the 
Lemaitre-Tolman (LT) metric to observations. One ap- 
proach uses low-z series expansions (TgI. [iTj. and a more 
general approach has shown how observational data fully 
determine an LT model [H, [H 113, HH, [H, |23, • A 
slightly different formulation can be found in [25[ , and an 
approach based on the characteristic initial value formu- 
lation of numerical relativity in [26j . All approaches must 
ultimately be implemented as numerical procedures. 

In this paper we provide a complementary approach to 
the observational metric. In the observational cosmology 
papers, the approach to solving the problem has focussed 
on using the observational data to analytically determine 
the metric functions. Here, with an eventual numerical 
scheme in mind, our approach emphasises firstly the full 
formal solution of the field equations for the observer 
metric, particularly noting the 4 arbitrary functions that 
emerge in the process and how the evolution is deter- 
mined by them, and secondly the algorithm for determin- 
ing the metric from observational data, especially show- 
ing how the arbitrary functions are fixed by the data. 
The first provides a better understanding of the geome- 
try and dynamics of the model, and the second shows the 
relationship between the data and the particular charac- 
teristics of the solution metric. We also give the explicit 
transformation between the OC and LT forms of this 
metric. 

Although spherical symmetry about the observer is a 
strong assumption, we regard it as a first step — a useful 
and important one — towards the more general case. 
A proper understanding of this simpler case is essential 
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for working with the more general forms of the observer 
metric. See [1^ [13] for a more detailed justification. 



II. SPHERICALLY SYMMETRIC OBSERVER 
METRIC 

We choose coordinates = {w, y, 9, 0), and we assume 
(a) spherical symmetry about the origin, (b) the observer 
is at the origin, and (c) the {6, 4>) surfaces are orthogonal 
to the {w, y) surfaces. We work in geometric units. With 
these, the metric is 

ds^ = -A^ dv? + 2AB dw dy + dVi^ , (1) 
where dV? = dO^ + sin^O dcf^ , (2) 
A = A{w,y), B^B{w,y), C = C{w,y). (3) 

Here C is an areal radius, and the lack of a dy^ term 
ensures that the constant w surfaces are null. 



dw — dfl 



ds^ ^ 



(4) 



We further assume that the matter is a zero-pressure per- 
fect fluid, comoving with the y coordinate. 



where the constant y, 
vectors such that 



(5) 

^, curves have timelike tangent 

{-A,B,0,0), u^ua^-l. 

(6) 



Here p is the proper density relative to observers on the 
comoving worldlines, m". 

We note that the null tangent vector ka within the 
constant w surfaces. 



fc"fca=0, (7) 



must be geodesic since it represents radial light rays, 
k'^Vak'' = 0, and in fact, for any K{y), ka = S^xly) 
is geodesic. Similarly, the dust particles should follow 
geodesies, and u°'V = leads to the two equations 
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"d,.,u^ ^^[1^ + ^ + 1^] [u-'f (8) 

(9) 

where A^ — dA/dw, Ay — dA/dy, etc. The second of 
these imposes a restriction of the metric functions, 

SU, = -Ay , (10) 

which reduces the first to {dwU^)/^^ — ~{dwA) / A, in 
agreement with the normalisation condition ([5]). 

Near the central worldline there is a spherical origin, 
where C — *■ 0. The origin conditions for this metric, 
giving the limiting behaviours of A, B and C near an 
origin, have been presented in several of the observational 
cosmology papers. 



A. Solving the EFEs 

The Einstein field equations (EFEs) G"^ = kT"'' - Ag"'' 
for this metric are 
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where k = Stt, and the conservation equations, VfcT"'' 
are 

V7 rrwb Pw , P f^Cw 25„, Ay 



A^ A^ \ C B B 



^ =0 (15) 



^ rj^yb _ P{B^ + Ay) _ 



(16) 



Not surprisingly, we can obtain (|10p directly from ([T 
since we don't expect the density to be zero. 
From & above we obtain 
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where two terms cancelled because of (|10p . This can be 
written as 



d ( C,„ C,, 



dw\ A B 
which solves to give 



= 



(18) 



(19) 



where W(y) is an undetermined function of integration. 
Next, from (fT9|) we have 

Cy=B{w-^ 
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Cwy — I W j i^ACww AyjCw) , (20) 



which combine with the G^^ equation (fT3|) to give 



-C^iW^ ^1) = C^C^A , (21) 
where ([TU]) was used again. The solution here is 



CC^ 



(22) 



^2 



C{W^ ~ 1) 



2M(y) , (23) 



where M{y) is a second undetermined function of inte- 
gration. 

Equations (|lip and ^3]) give the same as the y deriva- 
tive of (1231): 



KpBW = (G""" -I- A5""")^^BW - (G'"^ Ag'"^)^BG 
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where and its y derivative were used. Therefore the 
density is given by 



Kp 



2My 

CPBW ' 



(25) 



which clearly satisfies (115]) . Similarly, the Kretschmann 
scalar is 



A — it Uabcd 
12M?, 
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(26) 



The solution (|23p can be re-written as an evolution 
equation for G: 



G„ , 2M , , AG2 



(27) 



where f{y) = W^ -1 ^ W ^ \/l + f , 

(28) 

and the sign depends on whether G is increasing or de- 
creasing with time. In addition, equations (fT9|) and ([27|) 
give 



Eq (|27|) is clearly allied to the LT evolution equation 
(|A2[) . except that it contains two unknown functions, G 
and A, so it cannot be solved as is. In this paper, where 
a function is transformed between coordinates, we write 
e.g. C{'w,y) = C(t,r), meaning the two forms have the 
same numerical value at any given event, but they don't 
have the same functional dependence on their arguments. 
To solve (P7)) , we define t along the worldlines of constant 
y by 



t 



Adw 



const y 
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which converts (|27l) to 



dC 



± 



+ / + 



dt = t — a{y) 



(30) 



(31) 



In principle this gives us t{C, y) or C{t, y), and introduces 
a(?/), the initial t value at each y, as a third free function 
of integration. Clearly t is the proper time along each 
worldline, and the solutions to ([5T|) are identically those 
of the LT metric. When A = 0, the solutions for each of 
the cases / > 0, / = and / < are well known, and are 
often given parametrically, {C{ri,y),t{r/,y)}. However, 
we don't yet have a transformation between t and w, 
since by ((30)) 



t = / Adw 

' const y 



A tn, 



(32) 



we must know A to calculate w and vice versa. Now from 
(fTU| we find 



B„ 



Ay — iwy 



B = ~ty + P{y) . (33) 



The function /3(y) reflects a freedom in the definition ((30)) 
oft, 

t -> t + a{y) , 13 ^ P-ay , (34) 
which we shall remove. We next define 

r = r{w, y) = y , = , ry = 1 , (35) 

and in the next few equations we use r when it is paired 
with t, as in Cr = drC{t, r), but y when paired with w, 
as in Cy = dyC{w, y). By requiring that our t coordinate 
be orthogonal^ to r, viz: 



= 5°' [dat) [dbr) = 
A + Q iP~B 



^y ^y '^w ^y ^y 



AB 



AB B^ 

^ /3 = , (36) 



AG2 



(29) ^ Without this, /3 and several subsequent equations contain 
large extra terms. 
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we reduce the freedom in t to a constant translation, i.e. 
a is a constant, but henceforth we shaU drop it from our 
equations. Evidently, then, B is the negative of the rate 
of variation of proper time with respect to y down the 
past null cone. Even though wc know t{C,y), we can't 
calculate ty unless we know how to hold w constant. Now 
the transformation between C{w, y) and C{t, r) allows us 
to write 



Cy City + CrTy 



-CtB + Cr 



(37) 



where Ct = dC/dt and Cr = dC/dr. Combining §7\ 
with 111) leads to 



B 



VTTf 



(38) 



and since C{t,r) is known, this gives us B{t,r). Using 
|), ([M]) and we obtain the differential equation 



VTTf 



(39) 



This equation specifies how much t changes for a given 
y change, when w is constant, so it may be integrated 
down the null cones, i.e. along constant w, from the ori- 
gin outwards, giving t(w,y). The boundary conditions, 
fixed say at the origin y = 0, give us a 4th undetermined 
function, 7(w) = t{w,0). Actually, this fixes the varia- 
tion of w with respect to t, as t is fixed by integrating 
([3T|l . An obvious choice is w — t\o — j- Having solved 
we can then convert C{t,r) to C{w,y) using 



C{w,y) = C{t{w,y),y) , 
(40) 



C{t,r) plus t{w,y) 



and we finally determine A{'w, y) from (|32p . and B{w, y) 
from (|33p , or possibly (|38p . The algorithm for calculating 
the model evolution is detailed in section HVFI 

Having completed the solution, we see that this metric 
has 4 arbitrary functions — f{y), Mijj), a{y) & 7(w) — 
of which 7(w) represents a freedom to rescale w that is 
most naturally set to 7 = w. The physical meanings of /, 
M and a are exactly as in the LT model; they represent 
two physical relationships plus a freedom to rescale y. 



III. OBSERVABLE QUANTITIES 



Let wq label the past null cone (PNC) of present day 
observations by the central observer at {t,r) = (to,0) = 
(w,y) — {wo,0), and let the evaluation of any quantity 
Q{w,y) on this PNC be denoted Q — [QJa = Q{w{),y). 
We assume that emitters follow comoving worldlines j/e, 
and the observer is at the central worldline, yo — Q- Let 
the evaluation of a quantity at the observer and the emit- 
ter be denoted Qo{w) = Q{'w,Q) and Qe{w) = Q{w,ye) 
respectively. However we will often drop the subscript e. 

The redshift of comoving sources on that null cone is 
given by the ratio of the light oscillation periods T mea- 
sured at the observer, o, and the emitter, e. 



(1 + ^) 



Ao dw 
Af, dw 



A 



A,. 



(1 + ^) 



(41) 



and we can put A^ = 1 since A^, — A{'Wq, 0) = dw^\wo = 
1 is the natural choice. 

The diameter distance of a source is the true diameter 
D divided by the angular diameter ^, and in a spherical 
metric it corresponds to the areal radius evaluated on the 
PNC, i.e. 



— djj — Ce 



(42) 



Similarly, if the absolute luminosity of a source is L, the 
apparent luminosity is i (or m and in the apparent and 
absolute magnitude), and dio is 10 parsecs, then the lu- 
minosity distance is 



^)/5 d 



10 



(43) 



By the reciprocity theorem [27|, l28|, |29|, may be con- 
verted to do using z. 



do ^dL{l + zf 



(44) 



In redshift space, {z,6,(j)), let n{z) be the density of 
sources, that is the number per steradian per unit redshift 
interval^. Suppose that there are dN sources in solid 
angle diu — sin 9 dO d(j) between redshift z and z + dz, and 
that /i(z) is the mean mass per source'^, then the mass in 
that volume element of redshift space is 



dAi = /i dN = iindw dz 



(45) 



The primary observables for cosmological sources are 
redshift z, angular diameter S, apparent luminosity £, 
and number density in redshift space n. Associated with 
each of S, £, and n is a source property, true diameter 
D, absolute luminosity L, and mass per source fj,, that 
is needed to relate the observations to the theoretical 
model. The observables all depend on redshift z, and 
since sources evolve with time, so do the source proper- 
ties. We assume that the source properties are known 
from a combination of observations and source evolution 
theories. 



The proper 3- volume enclosing these sources at the time 
of emission, as measured by comoving observers is 
spanned by 



dxl = 6^ dy 



dx^ 



S^d0 , 



dx-^=S-^d^, (46) 



^ Thus this n is different from the n used in the OC programme, 
which is number density on a constant time sUce, see section Ivl 
^ For a treatment with a variety of source types, see |20|| . 
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and evaluates to 
(fiv = rjabcd dx'l da^s — \/\g\ £0123 dxl dx\ dx\ 



(47) 



- A^B'^C^ sm^e\ ^dyded<t) 
= BC^ sin 6 de dcf) dy = BC^ duo dy , 

so that the mass in this fluid element is 

dM = pBC^duidy . (48) 

Clearly we have the following relationship between n and 



2 dy 



^in = pBC^ 

az 



(49) 



The apparent horizon is where C is maximum on any 
given constant w cone, 



Cy = 



(50) 



Now if the metric ([T]) is to be regular, and the density 
(P5)l and Kretschmann scalar (pS)) finite at such a point, 
then B must be non-zero. This would need the upper sign 
in — that is, the local matter-shells are expanding, 
C > — and 



2M , AC2 , 

6M + AC^ - 3C = , 



(51) 



along this locus. When A = this simplifies to C = 
2M. As shown in [2l|, this locus has considerable 
observational significance. See also [isj . 



IV. DETERMINING THE SOLUTION FROM 
OBSERVATIONAL DATA 



Given the above observational data on the past null 
cone, that is A{z), C{z), and fj,n(z), the solution process 
must determine the arbitrary functions /, M, a and 7. 
Knowing these enables all spacetime quantities to be cal- 
culated and evolved, via the results of section III Al We 
envisage a numerical solution, so the emphasis here is on 
laying out a solution algorithm, rather than on formal 
integrals and functional dependence. 



A. Gauge choices 

The observational data must determine the physical 
properties of the model, but cannot restrict the coor- 
dinate freedoms. Therefore we will have to make some 
gauge choices in order to effect the solution. Firstly, we 
set w — t along the central worldline, which implies 



7(w) 



w 



A{w,0) = tni{w,0) = 1 



Ao = 1 



A^ 



(1 + ^) 



(52) 



Secondly, we need to set the freedom in the y coordinate. 
We consider two options below. The 'LT' option specifies 



ty = -l 



B = l, 



(53) 



as in many LT approaches. The OC papers choose 
A{wo,y) — B{wo,y) on the PNC, so in the 'OC option 
we choose 



B = A 



1 



(1 + ^) 



-1 



(1 + ^) 



(54) 



B. DE for y{z) 

The coordinate y is of course not observable, but we 
have to determine it first. We define 



dy 



(55) 



and along the PNC we re-write our equations in terms of 
z derivatives rather than y derivatives; for any quantity 
Q{w,y), 



0. 



Qy — and Qyy — 2 ^ — 3 



(56) 



Evaluating ([TTl) on the PNC, and using (gHl), we find 

AzCz , BzCz Czz , C^(p;, Kp.nBip 

-\ — — -\ = , (57) 

AC BC C Cip 2C2 

which, upon substituting for A^/A from (|¥T|) . leads to 
/ 1 Bz , Czz , KpnB(p\ 

(fiz = f \ -r, ^ ^ + — 1 — — I (58) 

where all derivatives are now total derivatives along the 
PNC. 

At this point we must fix the gauge in order to freeze 
out the coordinate freedom. The two options given above 
each convert (|55|) to an ODE for pivf) completely in terms 
of observables: 



OC: 



Czz 



Kfinifi 



(1 + 2) 4 2(1 + z)CC, 



(59) 

LT: (^2). = ^2 ( — ^ + ^ + . (60) 

Integrating (|59p or ([50]) followed by ([55]) yields fi{z) and 

yi{z) = I tp^{z) dz . (61) 
Jo 

This allows us to convert between functions of z and func- 
tions of y on the PNC. 
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C. DE for M{z) & W{z) 



From ^ on the PNC and we obtain 
where W is found by putting on the PNC, 



The two gauge choices give 



OC: 



W = 



2C,{l + z) 



1 - 



LT: 



iP2 

2C, 



2M 



2M 



AC2 



AC2 



(62) 



(63) 



2¥'2 



(64) 
(65) 



Together and ([53]) or ([55]) constitute an ODE for 
M{z) that also generates W{z). Note that ([M)) requires 
ipi from (|59p and (j65[) requires ip2 from (j60p . TechnicaUy, 
this is a first order Unear inhomogeneous ODE for M, so 
the formal solution is well known. In practice, the two in- 
tegrals involved would both have to be done numerically, 
so it is less work to solve the ODE directly in parallel 
with and using say a Runge-Kutta method. 



D. Obtaining a{z) 

From dSSI and ^ on the PNC we get 

4 = -cpB , (66) 
which, in the OC & LT gauges, simplifies to the ODEs 



OC: t,= ^ 

(l + z) 

LT: 4 = -(p2 , 



(67) 
(68) 



thus giving the worldline proper time, i{z) or i{y), on the 
PNC. The appropriate ipi must be used in each equation. 
From dSI]) on the PNC, we write 



dC 



±V^ + /+T 



(69) 



where / is given in (|28p and r is the proper time from 
the bang to the PNC along the matter worldlines, and 
after performing the integral at each z, we calculate 



a{z) — t{z) — t{z) 



(70) 



We now have M, W = ^1 + / , and a, and we've 
also used up the freedom to rescale y by directly or indi- 
rectly fixing ty. The only undetermined function is 7(10), 
though we have already fixed Aq — tw(wo,0) = 1, and 
we should extend this to 



^ — w 



E. Evolving off the PNC 



(71) 



In principle, no amount of data on the PNC is suf- 
ficient to determine the future evolution of any part of 
the spacetime, because new information can arrive along 
succeeding incoming light rays. But since we have al- 
ready assumed a dust equation of state, in order to get 
the arbitrary functions, it is not unreasonable to expect 
the worldlines continue their dust evolution into the fu- 
ture. This same assumption is tacitly made when fitting 
a Robertson- Walker model to observational data. 

Away from the PNC, z is no longer a useful variable, 
and we should rather use y. Also, the gauge choices don't 
give B off the PNC, so gauge-specific equations such as 
([59|) or (l68|) are not applicable. One may determine the 
full evolution of C using the algorithm below, based on 
the solution of section III Al In addition, (fTTj) together 
with ([25| provides a cross-check on the calculated prop- 
agation of the metric components. 

The remaining question is whether, given that we have 
initial data for C, B Sz A on the PNC, there is a better 
way to evolve C than integrating over the entire {t, r) 
domain twice, first calculating C(t,r) and then finding 
t(w,y) and converting to C^w^y). Can we integrate di- 
rectly with respect to w, giving C{w, y) straight off? This 
would be especially important if there were detectable 
time evolution in cosmological observables. The key dif- 
ficulty is that we don't know any of A, B or C away from 
the PNC and the central worldline, and though we have 
direct evolution equations for B^, and C^,, there isn't one 
for Aw. 

Once the arbitrary functions W , M and a are known, 
the observational data plus the gauge choices give us all 
of A, U & C on an initial constant w null cone, wq. For 
clarity of argument, consider Euler integration. Evolu- 
tion equations for C and B follow from ((27|) and (fTO|. 



Ci+i =Ci + {C^)^ dw , C„ = AV , 



(72) 



B,+i = B., + {Bw)^ dw , B^, = -Ay , (73) 

but the difhculty is finding an evolution equation for A. 
For example, (|19p . in the form 



A^+i 



+ 1 



W - {Cy\ + l/B,+ l W - (Cy), + l/B,+ l ' 



(74) 
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does not help because Ai+i cancels out. The yy EFE 
contains A^, but as soon as we substitute C^^, = A^^jV + 
AVw, then vanishes from the equation. The best we 
can do is to put — —Ay and = AV and Cww — 
A-wV + AVu] in the wy EFE, to obtain an expression for 
Ay /A that is free of w derivatives; but this requires an 
integration along constant w, and is at least as much work 
as solving ([55)1 . Similarly the 96 EFE gives an expression 
for dUAy/A). 

Combining ^ and (HH) with ^ eliminates Ak B, 
but leads to a second order non-linear PDF for C that 
depends only on M, f (or W), and their y derivatives. 
This would also require a double numerical integration 
over the {w, y) space. We have not been able to cast 
it in a simpler form, and we don't regard it as a better 
alternative, so we don't write it out here. 

Evidently there is no direct integration off the PNC 
along the constant y worldlines, though it should be pos- 
sible to program the numerical integration as a single 
sweep across the spacetime. 



F. The Algorithm 

The procedure for obtaining the observer metric from 
observational data may be presented as a two-part algo- 
rithm, the first for obtaining the undetermined functions 
from the data, and the second for calculating the model 
evolution from the functions. The arbitrary functions M, 
^ = \/l + / a-nd a are obtained as follows: 

• Assume the following observational data for a large 
number of sources on the PNC: 

redshift z, 

apparent luminosity i and absolute luminosity L, 
(and/or angular diameter 5 and true diameter D) 
number density of sources in redshift space n 
and mass per source /i. 
From these calculate diameter distance C — do{z) 
using ([^^ and and redshift space mass 

density fin{z). 

• Make a gauge choice, as in i jlV Al which fixes B 
and ty. 

• Integrate down the PNC one of equations 
(|58ll /([59 |) /([60 | . as appropriate to the gauge choice, 
which gives (p{z), then integrate (p{z) as in (|6ip to 
produce y{z) z{y). 

• Integrate ([5^ with the appropriate choice of 
(IMll/lEl/dSSl) down the PNC to calculate M{z) 
and W{z) -> M{y) & W{y). 

• Integrate the relevant choice of ([55]) to 
give the time on the PNC t{z); integrate ([69]) 
along each constant y worldline, producing t(z) the 
proper time from bang to null cone; then calculate 
the bang time a{z) from ([7(1)1 i{y) & a{y). 



Having found the 3 arbitrary functions, the evolution of 
the model is determined as follows: 

• Use equation (j3ip and integrate up and down each 
matter worldline to evaluate t{C,r) — > C{t,r) ev- 
erywhere (r = y). Initial conditions are provided 
on the PNC by C{z{y)) and i{y). In practice, this 
could be done in the same step as the r integration 
above. 

• Choose the gauge function ^{w) — t{w, 0) to fix w 
all along the central worldline. 

• Knowing C{t,r), calculate Cr everywhere; and 
hence find B{t,r) everywhere from psp . 

• From each w on the central worldline, integrate 
equation ([55)) to trace the {t, y) locus of its past 
null cone, allocating each point the same w, thus 
obtaining w{t,y) 1(111, y). 

• Calculate C{w,t) = C{t{w,y),y) and B{w,t) = 
B(t{w,y),y) everywhere, as shown in (|40p . 

• Differentiate t{w, y) to find A(w, y) according to 
([32]). 

The above steps are written for clarity rather than nu- 
merical efficiency. In coding it, certain steps may be com- 
bined. As explained in [23, [2^, the neighbourhoods of 
the origin, the bang, parabolic worldlines, and the max- 
imum in do require special numerical treatment. 

V. RELATIONSHIP TO OTHER WORK 

This solution must obviously be a version of the LT 
metric, so it would be useful to see the transformation. 
We give this in appendix \K\ 

We here compare the present paper with earlier work, 
particularly noting any differences, and we give the con- 
version between the different notations that have been 
used 

The notation for coordinates, (w, y, 0, 0); and metric 
functions A, B, C is common to all the papers, as are 
those for the redshift z and the primary PNC wq. In the 
OC papers, the w and y partial derivatives are written 
C ^ dujC ^ , C ^ dyC ^ etc. 

Concerning [j] , we note there is a factor of k missing on 
the right of their equations (13), or it has been absorbed 
into the definitions of /i, fiQ and p. We will assume the 
latter. They use a completeness factor F — the fraction 
of sources that are actually counted. Whereas survey 
completeness is a significant concern, in our paper we 
have assumed it has been corrected for. (The effect of 
systematic errors was investigated in {24].) Combining 
our ([T^ and ([TB]) gives 

Pw ( 2C^ B^ \ 
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p{y) 

C^B 



2My 

'kW 



(75) 



which is (21a) in d]. As noted in eq (30) of [J does 
not hold, so cqs their (34)-(45) are incorrect. 

In Q , they effectively obtain orthogonality of the con- 
stant t and y surfaces in their (13) by comparing the 
matter flow lines in the LT and OC metrics. Their vari- 
able iV* (y) is never interpreted — it is Btt times the total 
gravitational mass M within shell y, divided by the mean 
galaxy mass. Since M contains both the integrated rest 
mass and the curvature, N^: is not proportional to the 
number of galaxies (except near the centre). Item 4 in 
the corrigenda is misleading; to get that solution, put 
u — F — V into -\- 2uv = 1 — mN^^/C , multiply through 
by C, rearrange, differentiate with respect to y, and use 
the derivative of (27), giving 



{2CvF)' -f mFN' = {Cv^ + C)' 



(76) 



which can be evaluated on the PNC and integrated to 
give F . Note the prime in their paper (and in this para- 
graph) is d/dy while w is held constant. This must be 
remembered when taking the prime derivatives on the 
right of their (29), where functions are expressed in terms 
of i, y, and parameter T . In fact, we recommend inte- 
grating the negative of the left hand side of (29), to get 
the null cone path t(y), and then determining T rather 
than T' . Thus the derivative of the parametric solution, 
which includes calculating F', is not needed. There seem 
to be some oddly placed factors of 47r: if their (17) and 
the equation below their (10) are correct, then their p 
is 47r times the density. Also their (17) and (27) imply 
their mN,, is our M, whereas their (25) implies it is our 
2M. We will assume the latter, as their equations in- 
volving p do not play a significant role. Otherwise that 
paper (with corrigenda) is basically correct, though a lit- 
tle circuitous, and the propagation of the solution off the 
observer's PNC in observer coordinates is not discussed. 

In ^TTI, p and /io also contain absorbed factors of k. 
The function of integration l{y) in their (47) corresponds 
to I3{y) in our p3|l. but, without the orthogonality con- 
dition, they haven't set it to zero. In steps 2 and 3 of 
their integration procedure, they argue that A{w, y) must 
have the same functional form as A(?«o,y), though they 
acknowledge that this does not mean simple replacement 
of wq by w. In their FLRW example, they use the central 
conditions plus homogeneity to obtain 1 — > w/wq va. the 
formula for A. Unfortunately, we cannot see how this can 
be implemented in general, away from the origin, when 
we aren't assuming homogeneity. 

The table below summarises the correspondence be- 
tween the different notations. 
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VI. CONCLUSIONS 

We have presented two methods of solving the spheri- 
cal observer metric; firstly a formal solution of the EFEs 
in terms of arbitrary functions, and secondly a procedure 
for determining the arbitrary functions and the metric 
evolution from observational data. Both are given in a 
clear, step by step manner. We have been precise about 
where the gauge choices are made, and have considered 
two possible sets. Our main aim was to be able to lay 
out a solution algorithm that could be coded for numer- 
ical implementation. This naturally divides into to two 
distinct stages: determining the arbitrary functions from 
given observational data, and determining the spacetime 
evolution from known arbitrary functions. These are 
given in ijlVFI The understanding thus gained will be 
useful for generalisations away from spherical symmetry. 

We emphasise that, if we are given observational data, 
then they fix the arbitrary functions, and we are only free 
to make the gauge choices that pin down the coordinate 
freedoms. On the other hand, if we choose all the arbi- 
trary functions, then the observational relations are al- 
ready fixed, and there's no room to fit observations. We 
find previous solution methods have not so clearly dis- 
tinguished the two, and can be hard to follow. We also 
find those methods were sometimes ambiguous about the 
functional dependence of their solution functions, espe- 
cially when extending the solution off the PNC, so that 
it was not always apparent how to proceed with calcula- 
tions. 

The OC papers all start from the fluid ray tetrad equa- 
tions, though a number of different solution methods and 
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notations are used. Our main contribution is in clarifying 
how the model evolution is to be calculated. Although 
the error in 0] was corrected in Q , the latter does not re- 
ally address model evolution, and the method suggested 
in |Tl[ may not be practicable. Also, unlike [1, Q 
our method does not need to refer to a known equivalent 
metric to restrict the solution 

In the formal solution, our approach has a lot in com- 
mon with previous approaches up to ([33)) . but thereafter, 
up to the final solution ([iO]) . it is new. In the solution 
from observations, our method is similar to others in sec- 
tions [IVB] and to some extent HV D I where y{z) and a{z) 
are found, but differs in sections IIV CI and IIV El where 
M{z) and W{z) are found and the evolution off the PNC 
is discussed. In particular, the demonstration that the 
evolution from one constant w null cone to the next nec- 
essarily involves an integration down the null cone at 
each step, is new. The orthogonality condition leading 
to /3 = is important for the solution from scratch, but 
up to now has only been obtained indirectly by compar- 
ing with the LT metric. 

The available data on galaxy observations only extends 
to a finite redshift z. This is well suited to the algorithms 
presented here, which involve integrations along the null 
cones outwards from the centre, and integrations along 
the constant y worldlines. The LT functions are fully de- 
termined within the range of reliable data. Even if quite 



large scale inhomogeneity is discovered, it is still possi- 
ble the cosmos approaches homogeneity on even larger 
scales. 

In [3| it was shown that any reasonable C{z) and 
Hn{z) derived from observations could be fitted by an 
LT model with zero A. While this is true for nearly all z, 
[2ll | pointed out that the maximum in the diameter dis- 
tance is an exception. The properties of this particular 
locus allow the mass M to be calculated independently of 
([62)1 . With perfect observations, this cross-check would 
allow a determination of A. In practice, with real obser- 
vationSjthere would be systematic errors, and as shown 
in [2^, [ij] this cross-check may instead be used to detect 
and correct for systematic errors. Thus, one would need 
more than the observations considered here to distinguish 
non-zero A from inhomogeneity. 
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APPENDIX A: THE TRANSFORMATION OF LT 
TO NULL-COMOVING COORDINATES 



The LT (Lemaitre- Tolman) metric is [30, |3l| 



t\2 



1 + ./ 



dr^ + i?2 dn^ 



(Al) 



where R = R{t,r) and R = dR/dr. It depends on 3 
arbitrary functions, / = f{r), M = M{r) and a = a{r). 
The matter is comoving and has zero pressure. The evo- 
lution equation is 

Ai?2 



where R — dR/dr, and the density is given by 

2M' 



Kp ■ 



R^R 



(A2) 



(A3) 



The arbitrary functions each have a physical meaning; / 
gives the deviation of the constant t 3-spaces from flat- 
ness, and also gives twice the energy per unit mass of the 
dust particles; M gives the gravitational mass within the 
comoving shells of constant r, and a gives the local time 
of the big bang on each constant r worldline. 
We propose the transformation 



J 



t = t{w,y) , r = y 
d{t, r) 



d{w,y) 



tyj ty 

1 



, (A4) 



which retains y as a comoving coordinate, so the metric 
becomes 



ds^ 



-{tw dw -f ty dyY 



l\2 



jR') 
1 + / 



(r^ dw -\- ry dyY 



-R^dn'' 

-t^ dw^ — 2twty dw dy - 



R^ dn^ 



2 (R 



1 + ./ 



(A5) 
dy^ 
(A6) 



We want w to be a coming null coordinate, i.e. dw = = 

de = 



must give ds = 0, which leads to 

-R 



9yy * ^ 

ds^ = -tl, dw^ -\-2t 



' VT+f 
R 



VT+f 



(A7) 

dw dy + i?2 dn'^ , 
(A8) 



where the sign choice is because, on the past null cone 
(PNC), t must decrease as r = y increases. Eq (|A7|) is a 
PDE for t{w, y), and its solution will introduce a function 
of w, 7(w). 

Applying this transformation to the LT evolution equa- 
tion (ESI) 



T'lij - — - fi, tfi 



(A9) 



we find the new evolution equation in the new coordi- 
nates. 



2M , Ai?2 
Rw = ±tw^ — + f + —- 

Similarly, by transforming Ry we obtain 

R.' 



Ry - R ty 



R = Ry 



,VTTf T 



2M 



3 



(AlO) 

R 
(All) 

(A12) 
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where (IA10[) was used. The metric now becomes 
ds^ = -tl dw^ + 2R^t^_^_ 



i?2 dn^ 



(A13) 



Up to this point t(Tju, y) contains an undetermined func- 
tion of w, obtained when integrating (|A7p . We now spec- 
ify that, at the origin, y = r = 0, w is the observer's 
proper time, w = t, i.e. 

t{w,Q) = w, w{t,0)=t, t^{w,0) = l, 

wtit,0) = l. (A14) 

and this wiU fix the function of integration introduced in 
solving (|A7p . 



Comparing ([T]) and (lAlSp . it is clear that 



A = U 



B -~ty - 



in agreement with (15^ and 

The matter tensor transforms to 



ub _ / p/tl 



0, 



d I P^W ptyjty 

P ^10 P^y I 



rp rp rC jd 



(A15) 



(A16) 



